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ABSTRACT 

We report on the Chandra X-ray Observatory ACIS-S3 imaging observation of the Galactic globular 
cluster M4 (NGC 6121). We detect 12 X-ray sources inside the core and 19 more within the cluster half- 
mass radius. The limiting luminosity of this observation is « lO'^^ erg s^^ for sources associated 
with the cluster, the deepest X-ray observation of a globular cluster to date. We identify 6 X-ray 
sources with known objects and use ROSAT observations to show that the brightest X-ray source is 
variable. Archival data from the Hubble Space Telescope allow us to identify optical counterparts to 
16 X-ray sources. Based on the X-ray and optical properties of the identifications and the information 
from the literature, we classify two (possibly three) sources as cataclysmic variables, one X-ray source 
as a millisecond pulsar and 12 sources as chromospherically active binaries. Comparison of M4 with 
47 Tuc and NGC 6397 suggests a scaling of the number of active binaries in these clusters with the 
cluster (core) mass. 

Subject headings: binaries: close — globular clusters: general — globular clusters: individual (M4, 
NGC 6121) — novae, cataclysmic variables, contact binaries — X-rays: stars 



1. INTRODUCTION 

Faint X-ray sources (Lx ~ 10^^ erg s~^) were first 
discovered in observations made with the Einstein and 
ROSAT observatories (Hertz & Grindlay 1983; Ver- 
bunt2001). It was only with Chandra, however, that a 
large number of faint sources was identified: more than 
one hundred in 47 Tuc, and up to a few dozen each in e.g. 
NGC 6397, NGC 6752, NGC 6440, and wCen (Grindlay 
et al. 2001a,b; Pooley et al. 2002a,b; Rutledge et al. 2002). 
These faint sources represent a mix of objects with differ- 
ent X-ray luminosities. Brightest are the neutron stars 
accreting at a low rate from a companion (quiescent 
low mass X-ray binaries or qLMXBs), followed by white 
dwarfs accreting from low-mass companions (cataclysmic 
variables or CVs). The active binaries (ABs) tend to 
be the faintest, while radio pulsars with short periods 
(millisecond or recycled pulsars, MSPs) have X-ray lu- 
minosities in a similar range as CVs. The three types 
of chromospherically or magnetically active binaries are 
detached binaries of two main sequence stars (BY Dra 
systems), detached binaries of a main sequence star and 
a giant or a sub-giant (RS CVn systems) and contact 
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binaries (W UMa systems). 

Globular clusters contain many more neutron star bi- 
naries per unit mass than the galactic disk. Hence, if 
in a globular cluster one finds a binary with a neutron 
star (which is rather difficult to make via ordinary evo- 
lution of an initial binary) , it is highly probable that this 
binary was formed via a close encounter between stars 
(Fabian et al. 1975; Hills 1976). Isolated millisecond pul- 
sars in globular clusters are thought to have been formed 
in such binaries, and thus are probably also a result of 
close stellar encounters. 

On the other hand, binaries that are common in the 
field, are more likely to be of primordial origin in glob- 
ular cluster, including the chromospherically active bi- 
naries. Cataclysmic variables are relatively common in 
the Galactic disk, and those in globular clusters could in 
principle originate from primordial binaries or stellar en- 
counters. Pooley et al. (2003) showed that the number of 
faint sources above the threshold > 4 x 10^° erg s~^ 
in a cluster scales with its collision number P, which is 
a theoretical estimate for the number of close encoun- 
ters. Since the majority of such sources are cataclysmic 
variables, this suggests that most cataclysmic variables 
in globular clusters are in fact formed via close stellar 
encounters. 

Most clusters investigated with Chandra so far have 
relatively high collision numbers (see e.g. Table 1 in Poo- 
ley et al. 2003); because the limit to which X-ray sources 
in these clusters could be detected is relatively high, most 
X-ray sources known in these clusters are either neutron 
stars (accreting or radio pulsars), or cataclysmic vari- 
ables. So far, a sizable number of X-ray sources has 
been identified as chromospherically active binaries only 
in 47Tuc (Edmonds et al. 2003a, b). In this paper, we 
discuss a globular cluster with a relatively low collision 
number, M4 (NGC 6121). This cluster is a relatively 
nearby cluster with a moderate absorption (d = 1.73 kpc, 
Av = 1.32, Richer et al. 1997). The core and half-mass 
radii of M4 are 49'.'8 and 3f65, respectively (Harris 1996). 
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We use these values throughout the paper. 

A priori, we would thus expect that the majority of X- 
ray sources in this cluster to be chromospherically active 
binaries. That chromospherically active binaries exist in 
this cluster is evident from optical studies. For exam- 
ple, Kaluzny et al. (1997) discovered a number of optical 
variables, including several contact binaries. 

The presence of a recycled radio pulsar PSR B1620— 26 
(Lyne et al. 1988) in M4 is remarkable, but can be ex- 
plained by noting that the small collision numbers of 
many globular clusters still add up, so that at least some 
of the clusters with small F should contain a binary with 
a neutron star (Verbunt & Hut 1987). PSR B1620-26 is 
in a 191 day orbit around a white dwarf of ~ 0.3 Mq; 
the binary is accompanied by a third object of planetary 
mass in an orbit of 100 years (Thorsett et al. 1999; 
Sigurdsson et al. 2003). 

2. X-RAY OBSERVATIONS AND ANALYSIS 

M4 was observed for 25.8 ks on 2000 June 30 with 
the Advanced CCD Imaging Spectrometer (ACIS) on the 
Chandra X-ray Observatory with the telescope aimpoint 
on the back-side illuminated S3 chip. The data were 
taken in timed-exposure mode with the standard inte- 
gration time of 3.24 s per frame and telemetered to the 
ground in faint mode. 

Data reduction was performed using the CIAO 2.3 soft- 
ware provided by the Chandra X-ray Center (CXC)^. We 
reprocessed the data using the CALDB 2.21 set of cali- 
bration files (gain maps, quantum efficiency, quantum ef- 
ficiency uniformity, effective area) without including the 
pixel randomization that is added during standard pro- 
cessing. This method slightly improves the point spread 
function. We filtered the data using the standard ASCA 
grades, and we excluded both bad pixels and software- 
flagged cosmic ray events. Intervals of background flaring 
were searched for, but none were found, hence we simply 
applied the good-time intervals supplied by the CXC. 

2.1. Source Detection 

The CIAO wavelet-based wavdetect tool was employed 
for source detection in both the 0.5-6.0 keV band and 
0.3-10.0 keV band. We detected 41 sources on the entire 
S3 chip in the 0.5-6.0 keV band; for the wavdetect set- 
tings we used, approximately one of those detections may 
be spurious. Of these sources, 30 lie within the half-mass 
radius of the cluster. 

We also searched part of the adjacent S4 CCD since 
part of the half-mass region fell on this chip, but no 
sources were detected in this area. One additional source 
was detected on S3 in the broader 0.3-10.0 keV band with 
three counts, two of which were between 8 and 10 keV. 
The significance of this source, as well as its possible 
membership of the cluster, is questionable. For complete- 
ness, however, we leave it in our source list (Table ^ as 
the last source. The sources are numbered according to 
detected counts in the 0.5-6.0 keV band. 

From the logA^ — log 5* relationships of Giacconi et 
al. (2001), we expect between 5 and 6 background sources 
within the cluster half-mass radius. However, as M4 is 
located towards the bulge of the galaxy {I — 351.0°, 
b = 16.0°, Harris 1996) this relation would underestimate 
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Fig. 1. — X-ray image of a 5' X 5' region of M4. Shown are 
all counts in the 0.5-6.0 keV range. The cluster center, core and 
half-mass radii are indicated by a cross and two circles, respec- 
tively (Harris 1996). The detected X-ray sources are marked and 
numbered, omitting the "CX" prefix. X-ray sources classified as 
cataclysmic variables are indicated with triangles, chromospheri- 
cally active binaries with stars, millisecond pulsars with squares 
and unclassified or ambiguous sources are indicated with circles. 
See Section |^ for the source classification. Also shown are three 
archival pointings of the HST /WFPC2 used to locate optical coun- 
terparts to the Chandra X-ray sources. 



the galactic contribution of background sources and pro- 
vide a lower limit only. We have analyzed ACIS-S3 data 
from the October 15, 2002 Chandra observation of the 
low-mass X-ray binary (LMXB) MXB 1659-298, which 
has I = 353.8° and b = 7.3°. In the 27 ks of this obser- 
vation the 1/4 sub-array of the S3 CCD shows 5 sources, 
other than the LMXB itself, in the 0.5-6.0 keV band. 
Scaling this up to the whole chip and taking into account 
Poissonian errors would give rise to some 20 sources 
within an area of the size of the M4 half-mass radius. 
Given that MXB 1659-298 is located about 10° closer 
to the galactic center this would overestimate the number 
of background sources towards M4. Though both limit 
determinations are very different in nature, we expect 6 
to 20 background sources within the half-mass radius of 
M4, and zero or one within the core. 

All sources are consistent with being point sources. 
However, CX5 and CX9 overlap somewhat. They were 
identified by wavdetect as two separate sources. It is pos- 
sible that they are actually a single extended object, but 
the double-peaked nature of the image lends support to 
the interpretation as two point sources. 

There is some evidence for excess emission on the ap- 
proximately 32 200 pixels inside the core radius. The 
residual number of counts (i.e., the total counts minus 
those due to the 12 point sources inside the core) in the 
0.5-6.0 keV band is 438 ± 20 counts, while we would 
expect 287 ± 15 background counts on the basis of the 
number of counts that we measure in a source-free region 



X-ray Sources in M4 
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outside the core. This excess emission could be due to 
either unresolved point sources or diffuse emission. 

A 5' X 5' portion of the Chandra exposure of M4 is 
shown in Figure ^ The detected sources are encircled 
and numbered. The cross indicates the cluster center, 
while the circles denote the core and half-mass radius. 



2.2. Count Rates 

We extracted source counts in the following bands: 
0.5-1.5 keV (X,oft), 0.5-4.5 keV (X„iod), and 1.5-6.0 keV 
(Xhard). The detected count rate was corrected for back- 
ground, exposure variations, and foreground photoelec- 
tric absorption. We make these corrections in order to 
produce an X-ray color-magnitude diagram (CMD) that 
can be compared to the X-ray CMDs that have resulted 
from Chandra observations of other globular clusters. In 
addition, however, attention must be paid to differences 
in detector responses and, of course, exposure times and 
distances. 

The background count rate in each band was estimated 
from a source-free region on the S3 chip outside the core 
and to the northwest. The density of background counts 
in each band (for the 25.8 ksec observation) is 0.0043 
counts pixel"^ in Xgoft, 0.0077 counts pixel"^ in Xmod, 
and 0.0047 counts pixel"^ in Xhard- The background 
count rate in the core may be somewhat higher, but even 
factors of a few greater than this estimate have negligible 
effects on our analysis. 

In general, the exposure variations among sources were 
at the '^7% level or less, but CX30 had an exposure which 
was 20% less than the others. To account for these vari- 
ations in exposure, we applied multiplicative corrections 
based on the ratio of the average effective area of the 
detector at the location of a source in each of the three 
bands to that in the same band of CX6, which had the 
highest average exposure. The individual effective area 
curves for the sources were made using the CIAO tool 
mkarf. The average effective area of the detector at the 
location of CX6 in each of the bands was 538 cm^ (Xsoft ) , 
451 cm2 (Xmed), and 376 cm^ (Xhard)- 

While the previous corrections were relatively minor 
(at the few percent level or less) , the correction for pho- 
toelectric absorption is appreciable for M4. The conver- 
sion of optical extinction to column density (Predehl & 
Schmitt 1995) gives a value of Nh = 2.36 x lO^^ cm'^. 
We investigated the effects of such an absorption on three 
characteristic spectra: a 3 keV thermal bremsstrahlung, 
a 0.3 keV blackbody plus power law with photon index 
of r = 2, and a power law with a photon index of F = 2. 
The effects were most prominent in the Xgoft band, where 
the absorbed count rate was a factor of 2.2-2.5 lower 
than the unabsorbed one (depending on the spectrum). 
Averaging the results of each spectrum in each band, 
we use the following correction factors: 2.38 (Xgoft), 1-87 
(Xmcd), and 1.10 (Xhard)- TablcQllists both the observed 
and fully corrected counts in each band. The effect of 
the absorption correction on the X-ray CMD (Fig. EJ is 
a uniform shift of all points 0.27 units on the left axis 
and 0.34 units on the bottom axis. The bottom and left 
axes give the X-ray color and magnitude without this 
shift (they do, however, include the small corrections for 
background subtraction and exposure variations) . 
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Fig. 2. — X-ray color-magnitude diagram. The X-ray color is de- 
fined as the logarithm of the ratio of Xgoft (0.5-1.5 keV) corrected 
counts to Xjjard (1-5—6.0 keV) corrected counts, and the magni- 
tude is the logarithm of Xmcd (0.5-4.5 keV) corrected counts. Our 
correction for photoelectric absorption has the effect of uniformly 
shifting the data +0.27 units on the left axis and -1-0.34 units on the 
bottom axis. The bottom and left axes provide the absorbed color 
and magnitude scales C'), i.e., the observed colors and magnitudes 
uncorrected for absorption. For the sake of clarity, only a few error 
bars (Icr) are shown. Sources CX14, CX18, CX28 and CX31 were 
not detected in either the Xgoft or Xji^rd band, respectively. We 
illustrate their limits on logXgoft/Xijard by adopting a single count 
in these bands. The X-ray sources are marked as in Fig. lU See 
Section|l]for the source classification. 



2.3. Spectral Fitting 

We used the CIAO tool dmextract to extract spectra of 
sources CX1-CX4 in the 0.3-10 keV range. We binned 
the spectra to have at least 10 counts (20 counts for CXI) 
per bin and fit them in XSPEC (Arnaud 1996) using x 
statistics. 

For CX2, CX3, and CX4, three different models (with 
absorption) were fit: thermal bremsstrahlung (TB), 
blackbody with a small (~20%) power-law contribution 
(BB+PL), and power law (PL). We fixed N/f to the value 
from optical extinction. As expected for such low-count 
spectra, very few fits could be formally ruled out. We es- 
timated the unabsorbed source luminosities by averaging 
the results from the three best-fit models for each source. 
The spread in Lx of the three models was ~15% for each 
source. Fitting a linear relation to these luminosities 
versus corrected X,„ed counts, we have estimated the un- 
absorbed luminosities for sources CX5-CX31 based on 
their Xmod counts (note though that any differences be- 
tween spectra are not accounted for). These are listed in 
Table El 

The X-ray color and magnitude of CXI (Fig. ^ sug- 
gested that it might be a CV, but the spectral fits using 
an absorbed TB model (with Njy allowed to vary) pre- 
ferred the maximum temperature allowed by the model 
in XSPEC, namely, kT = 200 keV, which gave a 
X^/d.o.f. of 36.0/39. Fixing the temperature at the more 
reasonable value of kT — 25 keV still allowed for a sta- 
tistically acceptable fit (xVd-O-f- of 42.1/40). The TB 
results indicated a rather hard spectrum, and an ab- 
sorbed PL model (again with N^f allowed to vary) gave 
a best-fit photon index of 0.99 ± 0.17 and best-fit N// of 
(2.6tj j) X 10^1 cm-2 with a xVd.o.f. of 33.3/39. The 
PL fit is shown in Fig. |2| and the unabsorbed luminosity 
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Table 1. M4 X-Ray Sources 



Detected Counts/Corrected Counts'^ 
R.A. Decl. Lx(0.5-2.5 keV) 



Source^ 




(J2000.0)'^ 






(J2000.0) 


b 




Xmcd 


-^hard 


(erg 


s 


,-l)d 


Counterpart^ 


IDf 


CX 1. . . 


16 


23 


34.128 


(1) 


-26 


31 34.85 


(2) 


154/381.2 


379/736.4 


275/313.0 


8.3 


X 


10^1 


Opt.?/R9? 


cv 


CX 2. . . 


16 


23 


46.399 


(3) 


-26 


31 15.67 


(8) 


48/117.5 


72/138.3 


30/33.8 


6.0 


X 


1030 


Opt. 


cv 


CX 3. . . 


16 


23 


38.073 


(3) 


-26 


31 38.18 


(4) 


46/113.6 


75/145.2 


29/32.8 


5.4 


X 


1030 


Opt. 


AB 


CX 4. . . 


16 


23 


34.321 


(4) 


-26 


30 39.27 


(5) 


28/68.8 


56/108.9 


36/41.3 


3.6 


X 


1030 


Opt. 


CV/AB 


CX 5. . . 


16 


23 


34.300 


(4) 


-26 


29 56.47 


(6) 


23/55.6 


35/66.7 


12/13.4 


2.5 


X 


1030 


V56 


amb. 


CX 6. . . 


16 


23 


38.058 


(9) 


-26 


29 21.81 


(14) 


15/35.2 


26/47.8 


12/12.9 


1.8 


X 




- 




CX 7. . . 


16 


23 


45.848 


(12) 


-26 


28 55.03 


(19) 


23/54.9 


25/46.3 


3/2.9 


1.8 


X 




Opt. 


fg- 


CX 8. . . 


16 


23 


31.478 


(9) 


-26 


30 57.84 


(7) 


17/41.2 


24/46.0 


7/7.8 


1.7 


X 




V52 


AB 


CX 9. . . 


16 


23 


34.224 


(5) 


-26 


29 55.69 


(5) 


17/41.0 


23/43.7 


6/6.7 


1.7 


X 




V56 


amb. 


CX 10. . . 


16 


23 


35.047 


(8) 


-26 


31 19.18 


(7) 


14/34.0 


23/43.9 


9/10.0 


1.7 


X 


103O 


Opt. 


AB 


CX 11. . . 


16 


23 


32.399 


(7) 


-26 


30 45.63 


(6) 


8/19.4 


18/34.7 


10/11.4 


1.3 


X 




- 




CX 12. . . 


16 


23 


38.205 


(4) 


-26 


31 54.21 


(6) 


16/39.3 


17/32.7 


2/2.2 


1.2 


X 


iU 


Rad./Opt. 


MSP 


CX 13. . . 


16 


23 


34.326 


(7) 


-26 


32 02.33 


(10) 


11/26.9 


12/22.9 


2/2.1 


8.7 


X 


1029 


Opt./V49 


AB 


CX 14. . . 


16 


23 


25.948 


(15) 


-26 


33 54.62 


(10) 


0/- 


12/23.1 


12/13.9 


8.7 


X 


1029 


- 




CX 15. . . 


16 


23 


36.769 


(8) 


-26 


31 44.60 


(11) 


10/24.6 


11/21.1 


1/1.0 


8.0 


X 


1029 


Opt./V48 


AB 


CX 16. . . 


16 


23 


33.682 


(6) 


-26 


34 17.27 


(10) 


2/4.7 


6/11.4 


7/7.9 


4.3 


X 


1029 


- 




CX 17. . . 


16 


23 


35.975 


(11) 


-26 


31 01.89 


(8) 


7/16.9 


9/17.0 


2/2.2 


6.4 


X 


1029 


- 




CX 18. . . 


16 


23 


45.767 


(9) 


-26 


31 16.85 


(14) 


7/17.0 


7/13.2 


0/- 


5.0 


X 


1029 


Opt./V55 


AB 


CX 19. . . 


16 


23 


28.953 


(8) 


-26 


29 51.47 


(9) 


4/9.5 


7/13.3 


3/3.3 


5.0 


X 


1029 


- 




CX 20. . . 


16 


23 


36.881 


(11) 


-26 


31 39.43 


(15) 


5/12.2 


7/13.3 


2/2.2 


5.0 


X 


1029 


Opt. 


AB 


CX 21. . . 


16 


23 


34.610 


(5) 


-26 


32 04.60 


(19) 


4/9.7 


7/13.3 


3/3.3 


5.0 


X 


1029 


Opt. 




CX 22. . . 


16 


23 


33.362 


(9) 


-26 


31 45.65 


(15) 


5/12.2 


6/11.4 


1/1.1 


4.3 


X 


1029 


Opt. 


AB 


CX 23. . . 


16 


23 


40.189 


(20) 


-26 


29 26.10 


(25) 


5/11.6 


6/10.8 


1/0.9 


4.1 


X 


1029 






CX 24. . . 


16 


23 


42.096 


(12) 


-26 


31 37.03 


(13) 


3/7.2 


5/9.4 


2/2.2 


3.6 


X 


1029 


Opt. 


amb. 


CX 25. . . 


16 


23 


33.503 


(9) 


-26 


32 30.28 


(8) 


4/10.0 


5/9.8 


1/1.1 


3.7 


X 


1029 


Opt. 


AB 


CX 26. . . 


16 


23 


38.884 


(11) 


-26 


31 48.26 


(8) 


3/7.3 


5/9.5 


2/2.2 


3.6 


X 


1029 


Opt. 


AB 


CX 27. . . 


16 


23 


33.290 


(5) 


-26 


31 57.81 


(25) 


4/9.7 


5/9.4 


1/1.0 


3.6 


X 


1029 


Opt. 


AB 


CX 28. . . 


16 


23 


34.969 


(12) 


-26 


32 24.65 


(10) 


4/9.9 


4/7.6 


0/- 


2.9 


X 


1029 


Opt. 


AB 


CX 29. . . 


16 


23 


19.486 


(19) 


-26 


31 43.39 


(26) 


1/2.3 


1/1.7 


2/2.3 


6.5 


X 


1028 






CX 30. . . 


16 


23 


28.396 


(9) 


-26 


30 22.25 


(23) 


2/5.9 


3/7.0 


1/1.4 


2.6 


X 


1029 






CX 31. . . 


16 


23 


29.424 


(9) 


-26 


30 00.19 


(11) 


1/2.3 


1/1.7 


0/- 


6.5 


X 


1028 







''Sources are numbered according to their total counts. 

'^The Chandra positions have been corrected for a shift of 0'.'15 in decHnation (see il3.2i . Positional uncertainties are given in parentheses 

and refer to the last quoted digit an d ar e the centroiding uncertainties given by wavdetect. 
•^Corrections are described in Sect. IT2I The X-ray bands are 0.5—1.5 keV (Xgoft), 0.5-4.5 keV (X^cd) '^'nd 1.5-6.0 keV (Xji^rd)- 
"^For sources CX1-CX4, Lx comes from an average of the best-fit models for each source (Sect.|^^. A linear relation between Lx and Xj^^d 

counts for these sources was derived and used to estimate Lx for sources CX5-CX31 based on their X,„(,d counts. Typical uncertainties in 

Lx are ~ 15 %. 

"^Type of counterpart (optical and/or radio) found and associations (if any) with previously reported sources. The "V" numbers refer to 
optical variables by Kaluzny et al. (1997) and Mochejska et al. (2002). R9 is an X-ray source detected by ROSAT (Verbunt 2001). 

^Classification of the sources; see Sect.|3 The abbreviation MSP stands for millisecond pulsar, CV for cataclysmic variable and AB for 
active binary. Sources for which the classification is ambiguous are abbreviated with amb. and the foreground source CX7 is classified as 
fg.. See text for details. 



from this model is 6.8 x 10'^^ erg s ^ in the 0.5-2.5 keV 
band. 

3. OPTICAL OBSERVATIONS 

Three fields located inside the half-mass radius have 
been observed with the Wide Field Planetary Camera 
2 (WFPC2) aboard the Hubble Space Telescope (HST), 
and are shown in Fig. ^ Deep observations (dataset GO- 
5461) of three fields were obtained to determine the mass 
and luminosity functions of the main and white dwarf se- 
quence (Richer et al.l995, 1997) and two of these fields 
(see Fig. ^ gray outlines) fall within the cluster half- 
mass radius. These two fields were imaged for 11 800 s in 
F336W (hereafter U336), 15 000 s in F555W (V555) and 
5 500 s in F814W (/814) and were re-observed in /814, 
with exposure times of 5 680 s, (GO-8153) approximately 
5 years later. This allowed separation of field stars from 
the low luminosity cluster members through proper mo- 
tion properties (Bedin et al. 2001). Finally there exist 



shallow (249 s) V555 and (249 s) Isi4 observations (GO- 
6116) originally obtained to search for the optical coun- 
terpart of PSR B1620-26. These observations (Fig. [H 
black outline) cover nearly the entire region enclosed by 
the core radius. Of the 31 sources listed in Tabled 18 
coincide with the three HST /WFPC2 fields. 

This section outlines the data reduction, photometry 
and astrometry of the HST /WFPC2 images. 

3.1. Data Reduction and photometry 

The single HST images and the association products^ 
were obtained from the ESO archive. The single im- 
ages were already calibrated, including full bias subtrac- 
tion and flat-fielding. The association products consist 
of co-added, cosmic-ray cleaned images and association 
files, which contain the shifts between the single images. 
These shifts were used as input for the further reduction 
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Fig. 3. — Chandra spectrum of the brightest source, CXI. Dat- 
apoints are shown as crosses in the top panel, and an absorbed 
power-law model is shown as a solid line; residuals from this fit are 
shown in the bottom panel. Both the data and model have been 
divided by the instrument effective area for the purpose of plotting. 



of the data, using HSTphot 1.1 (Dolphin 2000a). This 
package is especially written for the reduction and pho- 
tometry of HST images. 

First, the HSTphot task mask was used on each image 
to mask bad pixels and image defects. This ensures that 
these pixels are not used in the further analysis. The task 
crmask was used to remove the cosmic ray hits, which 
arc widely present in the separate images. For each filter 
the images were compared against each other to remove 
the cosmic ray events. The HSTphot task hotpixels was 
then used to mask the known hot pixels in the WFPC2 
detectors. 

To obtain photometry for each individual image we 
used the photometry task hstphot of HSTphot 1.1. This 
task takes as input all available images in each filter and 
produces a master list of positions and magnitudes for 
each star found. The positions and magnitudes are deter- 
mined by fitting a model point spread function (PSF), to 
each star. After all stars are fitted it calculates aperture 
corrections and corrects for the charge transfer efficiency 
effect. These corrections, together with the zero-points 
are described in (Dolphin 2000b). 

3.2. Astrometry 

To search for optical counterparts to the Chandra X- 
ray sources wc aim to place both the X-ray and the 
optical frame onto the International Celestial Reference 
System (ICRS). We use this approach to improve the 
absolute pointing accuracy of Chandra and HST, 0'.'6 
and I'.'O (Icr) respectively (Aldcroft et al. 2000; Biretta 
et al.2000). 

To place the HST images onto the ICRS we use two 
intermediate steps; first, we align a ground-based im- 
age onto the ICRS using an astrometric catalog, then we 
align the HST images onto the ground-based image. To 
bring the X-ray frame onto the ICRS we use three optical 
identifications of X-ray sources. 

3.2.1. Astrometry of the HST images 

A 2 minute y-band image, taken on February 21, 2002 
with the Wide Field Imager (WFI) at the ESQ 2.2 me- 



ter telescope on La Silla, was retrieved from the ESO 
archive and used to calibrate the _ffS'T/WFPC2 images. 
The WFI has an array of 8 CCDs, each CCD having a 
8' X 16' field of view, giving a total of 33' x 34'. M4 
is roughly centered on one chip, and, to minimize the 
effects of geometric distortion, an 8' x 8' sub-image, con- 
taining the entire area within the half-mass radius, was 
extracted from this chip. 

We found 115 stars on this sub- image that matched 
entries in the USNO CCD Astrograph Catalog (UCACl, 
Zacharias et al. 2000). Of these stars, 91 were not satu- 
rated and appeared stellar and unblended and were used 
to compute an astrometric solution, fitting for zero-point 
position, scale and position angle. Five outliers, having 
residuals larger than 0'.'2 were rejected from the fit. The 
final solution has rms residuals of 0'.'05 in both coordi- 
nates. 

We calibrate the HST frames with the astrometric so- 
lution of the WFI image. First, the pixel positions of 
the stars in the HST datasets are placed on a single 
metaframe using the geometric distortion corrections and 
relative chip positions / orientations found by Anderson & 
King (2003). Next, using the astrometric solution in the 
FITS header of the HST images, nominal celestial coor- 
dinates are computed for some 300-400 of the brightest 
stars. These positions are then matched with stars on 
the WFI image and their centroids arc measured. Simi- 
lar selection criteria for saturation and blending are used 
as with the UCAC stars. An astrometric solution is cal- 
culated from the calibrated celestial coordinates and the 
metaframe positions of the HST stars. Outliers having 
residuals larger than three times the rms residual of the 
fit are removed, and a new solution is computed. This 
process is iterated until convergence. On average after 
convergence the astrometric solution contained some 200 
stars with rms residuals of O'.'OT in both right ascension 
and declination. 

The Icr uncertainty in the optical positions is 0'.'12, the 
quadratic sum of the positional uncertainty in the tie of 
the WFI image to the UCAC (O'.'O?) and the uncertainty 
in the transfer to the HST frames (O'.'OQ) and the un- 
certainty in the tie of the UCAC onto the ICRS (0'.'02, 
Assafin et al.2003). 

By placing the HST and WFI observations onto the 
ICRS with the UCAC we effectively set these images at 
the epoch of the UCAC observation, which is 1999.4. 
This implies that any average proper motion between 
the epochs of the HST and WFI observations has been 
removed. 

The difference in proper motion between cluster mem- 
bers and background stars (Bedin et al. 2003) possibly 
explains why our uncertainty in the calibration between 
HST and WFI is larger than what we have found for sim- 
ilar calibrations of data from other clusters (e.g. Bassa 
et al. 2003). 

3.2.2. Astrometry of the Chandra frame 

We use three optical identifications in the HST data of 
the Chandra, X-ray sources to place the X-ray frame onto 
the ICRS frame of the HST/WFPC2 images: the white 
dwarf companion of PSR B1620-26 and two W UMa 
variables. 

At the UCAC epoch, 1999.4, the position of 
PSR B1620-26 is aj2ooo = 16'^23™38?2147(5), 5j2ooo = 




"555 '814 "^336 '814 ''555 '814 

Fig. 4. — Color-magnitude diagrams (CMDs) of the H5'T/WFPC2 observations of M4. The first three panels show CMDs from the 
GO-6116 dataset for the PCI, WF3 and WF4 chips, respectively (no X-ray sources coincide with the WF2 chip). The fourth and fifth 
panel show two CMDs of the data in both GO-5461 fields. All stars inside the 95% confidence error circles are marked, while the candidate 
counterparts are numbered. The X-ray sources are marked as in Fig. The candidate counterpart to CX24 is plotted twice, as it was 
observed in both fields of the GO-5461 dataset and had varied in brightness. The candidate counterpart was brightest during the observation 
of the field closest to the core of M4. 



— 26°31'53'.'95(4), where we use the position and proper 
motion from Thorsett et al. (1999). According to our op- 
tical astrometry, this position is compatible with a white 
dwarf, at an offset of -0'.'04 ± 0'.'09 (-0!003 ± 0!006) in 
right ascension and 0'.'07 ± O'.'IO in dccUnation for the 
GO-5461 dataset. (The uncertainty in this offset is the 
quadratic sum of the uncertainty in the optical astrome- 
try and the uncertainty in the corrected pulsar position.) 
With our accurate position (0'.'12) for the white dwarf, 
we confirm the identification by Sigurdsson et al. (2003) 
and Richer et al. (2003) of this white dwarf as the com- 
panion of the pulsar. The position of X-ray source CX12 
coincides with that of PSR B1620-26. 

Five optical variables found by Kaluzny et al. (1997) 
and Mochejska et al. (2002), roughly coincide with Chan- 
dra X-ray sources, three of which are in the HST field(s) 
of view. Comparison of the Kaluzny et al. (1997) find- 
ing charts with our WFI and HST images shows that 
the variables are indeed excellent matches to the X-ray 
sources. The three variables coincident with the HST im- 
ages are all identified as W UMa binaries. V55/CX18 is a 
single star on the HST images, and we are confident that 



it is the star responsible for the optical variability. How- 
ever, the variables V48 and V49 both are blends of four 
and two stars respectively. In the case of CX13/V49 one 
of the stars is clearly above the main sequence, indicating 
that this star is the variable. However, for CX15/V48 all 
stars of the blend lie on the main-sequence, and none of 
these can be securely identified as the variable. 

We therefore use the optical positions of the white 
dwarf companion to PSR B1620— 26 and the variables 
V49 and V55 to compute the shift needed to place the X- 
ray positions of CX12, CX13 and CX18 onto the optical 
positions. We have a total of 6 measurements, as CX12 
is in the field of view of datasets GO-6116, GO-5461 and 
GO-8153, CX13 only in that of GO-6116 and CX18 in 
those of GO-5461 and GO-8153. The weighted average 
shift is O'.'OO ± 0'.'05 in right ascension and 0'.'15 ± 0'.'05 in 
declination, well within the la uncertainty of the abso- 
lute Chandra pointing accuracy. We apply this shift as 
a boresight correction to the X-ray positions. 
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Fig. 5. — 5" X 5" finding cliarts for tiie candidate optical counterparts. The finding charts are constructed from the coadded V555 
"association" images. For sources CX2, CX12, CX18 and CX24, images from the GO-5461 dataset were used, images from the GO-6116 
dataset are used for the other candidate counterparts. The 95% confidence uncertainties on the Chandra positions are overlaid on these 
charts, while the candidate counterparts are indicated with an arrow. No arrow is present for sources CX15 and CX21 as no candidate 
counterparts were found. The pixel scale varies between sources coincident with WF chips and the PCI chip (CX3 and CX26). The 
greyscale of these images is chosen as such to enhance the visibility of the candidate counterparts. 



3.3. Identification of counterparts in the HST images 

We search for optical/UV stars inside 95% confidence 
radii of the Chandra sources positions. The la uncer- 
tainty in a source position is the quadratic sum of the 
positional uncertainty for the X-ray source (Table^, the 
uncertainty in the optical astrometry (0'.'12) and the un- 
certainty in the X-ray boresight correction (0'.'07). The 
95% confidence radius is a factor [— 21og(l — 0.95)]^/^ = 
2.448 larger than the la uncertainty in the position. 

Each of the X-ray source error circles was visually 
checked for stars that were not found during the initial 
run of hstphot. Stars that were missed were added by 
hand and the photometry task was executed again. 

The resulting photometry lists were used to create a 
set of color- magnitude diagrams (CMDs), shown in Fig- 
ure^ The data from the GO-6116 dataset is separated 
into three panels, one for each _ff5'T/WFPC2 chip con- 



taining X-ray sources, as the cluster main sequence is 
displaced by about 0.1 magnitude in V555 — /814 between 
the individual chips. Differential reddening towards M4 
(Cudworth & Rees 1990) is the likely cause for this. The 
data from the two GO-5461 fields, shown in the last two 
panels, did not display these displacements in color and 
the data from all four chips is plotted. 

We used the hstphot information indicating the good- 
ness of fit for a star to exclude the diffraction spike 
artifacts from the star list. These artifacts are widely 
present, especially in the GO-5461 dataset, which had 
the longest exposures. For more information we refer to 
the HSTphot 1.1 manual^o. 

All stars within the 95% confidence radii of the Chan- 
dra source positions are marked in Figure 01 Stars with 
colors abnormal with respect to the cluster main se- 

|http: //Hww.noao . edu/staf f /dolphln/hstphot/1 
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Table 2. Optical/UV Counterparts to Chandra X-ray Sources 
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Optical positions and magnitudes of the 
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counterparts. The optical positional is given by the Chandra value 
plus offset (in true seconds of arc) in each coordinate. Magnitude 
uncertainties are given in parentheses and refer to the last quoted 
digit. Stars that were saturated in a pass band are denoted with 
"sat." and stars that are not detected in a band are denoted with 
"und." . 

''The HST GO dataset number from which the position and mag- 
nitudes are determined. For CX3, CX24 and CX26 the first entry 
for dataset G05461 and G08153 refers to the field nearest to the 
cluster center. 

'^This measurement is likely in error due to the proximity of the 
star to the edge of the chip. 



quence or giant branch are identified as candidate coun- 
terparts and are numbered in Fig.^ Positional and color 
information for each candidate counterpart is tabulated 
in Table 121 Finding charts are shown in Figure |31 

3.4. Identification of counterparts in the WFI image 

Several of Chandra X-ray sources coincide with stars 
on the WFI image. Two of these stars, identified by 
Mochejska et al. (2002) as variables V52 and V56 are 
coincident with X-ray sources CX8 and the possibly 
blended source CX5-9, respectively. The offset of the 
optical positions with respect to the Chandra positions 
in Table [2 is 0'.'06 ± 0'.'20, 0'.'22 ± 0'.'21 for CX8 and 
V52, -0'.'29 ± 0'.'22, -0'.'23 ± 0'.'21 for CX5 and V56 and 



0'.'73 ± 0'.'20, 0'.'55 ± 0'.'21 for CX9 and V56. 

The X-ray source CX7 coincides with the UCACl star 
22560073, a foreground star according to its proper mo- 
tion (Cudworth & Rees 1990, star A330 with V = 12.68). 
The offset of the corrected X-ray position, as listed in Ta- 
bleUI from the optical position from UCAC is 0'.'38±0'.'18 
in right ascension and — 0'.'16 ± 0'.'20 in declination. The 
error in this offset is dominated by the uncertainty in 
the X-ray position of CX7 in the Chandra frame; much 
smaller contributions to the error are due to the uncer- 
tainties in the X-ray to optical shift and in the optical 
position. Redetermining the X-ray to optical shift includ- 
ing CX7 with the other three optically identified X-ray 
sources, does not significantly change the value of the 
shift. 

4. SOURCE CLASSIFICATION 

In attempting to classify the X-ray sources detected 
in M4, we note that five of our optical counterparts to 
Chandra sources have a previously known classification. 
CX12 is a millisecond pulsar with a white dwarf com- 
panion (Thorsett et al. 1999); the optical counterpart 
corresponds to the white dwarf (Sigurdsson et al. 2003, 
see also Fig. p. CX13/V49, CX15/V48 and CX18/V55 
are contact binaries, and CX8/V52 is a BY Dra system, 
i.e. these four sources are all magnetically active binaries 
(Kaluzny et al. 1997; Mochejska et al. 2002). The optical 
counterparts to CX13 and CX15 as identified from the 
ground based data, have both been resolved into multiple 
objects with ir5'T/WFPC2. Kaluzny et al. (1997) noted 
that the optical counterpart to CX13, was brighter than 
expected for a contact binary with its period, placed at 
the distance and reddening of M4. This discrepancy is re- 
moved as our counterpart to CX13 is about a magnitude 
fainter than the blend (we note that the optical coun- 
terpart to CX13 appears to be too bright to be on the 
binary main sequence as it is about a magnitude brighter 
than a star on the main sequence with the same color. 
Figure 21). For the optical counterpart of CX15 the op- 
posite appears to be the case. It was seen as a cluster 
member by Kaluzny et al. (1997), but HST ob- 
servations indicate that V48 is a blend of 3 or 4 stars. 
Hence, the contact binary is too faint to be a member, 
and may be a background object. 

The possible blend of CX5 and CX9 coincides with 
the variable giant V56 found by Mochejska et al. (2002). 
If the giant is a member of a RS CVn binary then its 
variability and X-ray emission can be explained by mag- 
netic activity. However, without more information on 
the optical properties of V56 the classification remains 
ambiguous. 

In classifying the remaining X-ray sources we first look 
at the X-ray emission itself. A neutron star accreting at 
a low rate from a companion star, i.e. a low-luminosity 
(quiescent) low-mass X-ray binary (qLMXB), is charac- 
terized by a soft X-ray spectrum (black body color tem- 
perature < 0.3 keV) and a luminosity > lO'^^ erg s~^ 
(Verbunt et al. 1994; Rutledge et al. 1999). None of 
the X-ray sources in our sample shows this character- 
istic, and we conclude that M4 does not contain a low- 
luminosity low-mass X-ray binary with an accreting neu- 
tron star. 

We extract further information from the location of the 
optical star in the color-magnitude diagrams of Fig. 01 
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Fig. 6. — The X-ray luminosity Lx (0.5-2.5 keV) and the ab- 
solute V555 (approximately V-band) magnitude for X-ray sources 
with identified optical counterparts. Five types of X-ray sources 
are shown, quiescent LMXBs (diamonds), MSPs (squares), CVs 
(triangles), ABs (stars) and unclassified (crosses). The closed 
points are sources from 47Tuc (Grindlay et al. 200f, Edmonds 
et al. 2003a) and the open points are from NGC 6752 (Pooley et 
al. 2002). The numbered points are the optical counterparts to 
the Chandra X-ray sources in M4, the number corresponds to the 
source number. The absolute magnitude of the optical counterparts 
is computed from the V555 magnitudes and the V-band distance 
modulus as in Harris (1996) for 47Tuc and NGC 6752 and that of 
Richer et al. (1997) for M4. The dashed line of constant X-ray 
to optical flux ratio roughly separates cataclysmic variables from 
active binaries. 



stars bluer than the main-sequence are possible cata- 
clysmic variables, stars above the main-sequence possi- 
ble magnetically active binaries. The ratio of X-ray to 
optical flux is also useful, as it roughly separates the 
cataclysmic variables from magnetically active binaries 
(Verbunt & Johnston 2000; Pooley et al. 2002a). We il- 
lustrate this in Fig. IHl where we plot data from 47 Tuc 
(Edmonds et al. 2003a,b) and NGC 6752 (Pooley et al. 
2002a), and add our data from M4. 

CX2 is a probable cataclysmic variable, because it is 
blue (Fig. 01) and has a relatively high X-ray to optical 
flux ratio (Fig.|HJ). High X-ray to optical flux ratios sug- 
gest that CXI and CX4 also are cataclysmic variables, 
even though they are located on or close to the main 
sequence in the color- magnitude diagram (Fig. 0J . 

The X-ray luminosity of CXI is too high for a magneti- 
cally active binary of two main sequence stars (i.e. a BY 
Dra system), and its optical magnitude excludes that it 
is a magnetically active binary with a sub-giant (i.e. a 
RS CVn system). This, and the relative hardness of the 
X-ray spectrum (Fig.|3l indicate that it is a cataclysmic 
variable. We note that the optical counterpart to CXI 
is offset from the Chandra position by 2(7 and might be 
a chance coincidence. This would indicate that the ac- 
tual optical counterpart is even fainter, possibly lost in 
the glare of the bright star, leading to higher X-ray to 
optical flux ratios, as seen with most other CVs. 



The less luminous counterpart to CX4 is probably a 
cataclysmic variable; its location near the main-sequence 
in ^555 ~ -^814 has precedents in other globular clusters 
(Edmonds et al. 2003a,b). 

The optical counterpart of CX24 varies by more than 
a magnitude in brightness (Fig. 0)). It is located fur- 
ther above or to the right of the main-sequence than 
the binary sequence. Hence we believe that CX24 is ei- 
ther a fore or background object and unrelated to M4. 
The same argument could be made for CXIO. However, 
stars at similar positions in the color-magnitude diagram 
of other old clusters are X-ray sources and confirmed 
members (albeit unexplained, see e.g. Mathieu et al. 2003 
for M67 and Orosz & van Kerkwijk 2003 for NGC 6397) 
or probable members (Albrow et al. 2001, Edmonds et 
al. 2003a for 47 Tuc). We therefore consider CXIO a prob- 
able cluster member. 

Based on positional coincidence we identify CX7 with 
the UCACl star 22560073. This star has V = 12.68, B~ 
V = 0.88 (Cudworth & Reesl990), which is compatible 
with a K2V star at a distance of about 200 pc. Scaling 
the X-ray luminosity from Table from the distance of 
M4 to a distance of 200 pc gives w 2.4 x 10^^ erg s~\ 
well within the range observed for ordinary K2V stars 
(Verbunt 2001). 

The relatively low X-ray to optical flux ratios of the re- 
maining optical counterparts found with HST /WFPG2 
suggests that they are active binaries. The counterparts 
of CX13 and CX18, already known to be contact binaries 
from ground based data, are among these. (The V555 im- 
age of the optical counterpart of CX18 is over-exposed; 
we use the estimate V555 — (C/336 + /8i4)/2.) The oth- 
ers (CX3, CX20, CX22, CX25, CX26, CX27, CX28) are 
at or above the main-sequence, as expected for binaries, 
and we conclude that all of them are magnetically active 
binaries. The X-ray-to-optical flux ratio of any of the 2 
candidate counterparts of CX21 is comparable to that of 
an active binary. 

5. DISCUSSION 

The X-ray luminosities of the Chandra sources in M4 
are amongst the lowest ever observed in a globular clus- 
ter. More than half of the X-ray sources in this cluster 
have Lx < 10^° erg s~^, compared to 10% of the sources 
in 47 Tuc (Grindlay et al. 2001a) and 28% of the sources 
in NGC 6397 (Grindlay et al. 2001b). 

It is therefore not surprising that only one of the X-ray 
sources detected with Chandra may have been detected 
before. The position of the marginal ROSAT HRI source 
R9 (Verbunt 2001) is 7'.'5 from that of CXI. Even though 
this is further than expected from the error given for the 
ROSAT HRI source, we think the identification is prob- 
able. The X-ray luminosity of the ROSAT HRI source 
is only 1.3 x 10'^^ erg s~^, a factor six below the lumi- 
nosity detected with Chandra. The upper limit from the 
ROSAT PSPC observation is 1.5 x lO^^ erg s"^ (These 
luminosities are re-computed for a distance of 1.73 kpc 
and for the spectrum of CXI, and thus differ from those 
given by Verbunt; mainly because he used a distance of 
2.2 kpc.) This large variability indicates that CXI is a 
cataclysmic variable rather than a radio pulsar. Note 
that even if the ROSAT source is not identical to CXI, 
it still provides an upper limit to the flux of CXI during 
the ROSAT observations, and thus proves that CXI is 
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Table 3. Scaling Parameters of M4, NGC6397 and 47 Tuc 



Cluster 


log Po 


rc 


d 


Mv 


r 






(Lq pc-3) 


(") 


(kpc) 








M4 


4.01 


49.8 


1.73 


-6.9 


1.0 


1.0 


NGC 6397 


5.68 


3.0 


2.3 


-6.6 


2.1 


0.024 


47 Tuc 


4.81 


24.0 


4.5 


-9.4 


24.9 


12.4 



Note. — Values for central density (po), core-radius (rc), dis- 
tance id) and absolute visual magnitude {My) originate from Har- 
ris 1996 (version of February 2003). For M4, the values of po and 
My are computed for the distance and reddening of Richer et 
al. (1997). The collision number is computed from F oc pj'® 
and the core mass from Mc oc po r^. Values for F and Mc are 
normalized to the value of M4. 



highly variable. 

Several other X-ray sources coincide with optical vari- 
ables discovered by Kaluzny et al. (1997) and Mochejska 
et al. (2002). The sources CX13, CX15 and CX18 are 
coincident with W UMa variables V49, V48 and V55, 
with orbital periods of 0.283, 0.298 and 0.311 days, 
respectively. The upper limit on the X-ray luminos- 
ity of the other W UMa binaries found by Kaluzny et 
al. (1997), V44, V47, V50, V51, V53 and V54, is about 
3 X 10^^ erg s^^. Both detections and upper limits of 
these W UMa variables are in agreement with the range 
of X-ray luminosities of W UMa binaries in the ROSAT 
AU Sky Survey (St§pieh, Schmitt & Voges2001). 

CX8 coincides with the variable V52, which is classified 
as a BY Dra system with a period of 0.777 days. The 
possible blend of X-ray sources CX5 and CX9 coincides 
with the variable V56. No period is known for this object, 
but its variability (a rise of 0.1 mag in 4 days) and its 
location on the giant branch (Mochejska et al. 2002) are 
suggestive of an RS CVn variable. If it is, it is the first RS 
CVn binary detected in X-rays in a globular cluster - a 
marked contrast with old open clusters, where RS CVns 
dominate the X-ray sources (Belloni et al. 1998). The 
exposure of the eclipsing binary V54 is not as good as 
for the other sources, hence the upper limit on the X-ray 
luminosity of this source is somewhat higher. Detections 
and upper limits of these binaries are well within the 
range observed for magnetically active binaries, e.g. in 
the ROSAT All Sky Survey (Dempsey et al. 1993). 

In trying to determine the X-ray luminosity function 
of M4, we can use the number of detected sources as a 
function of luminosity. However, the number of excess 
counts, not allocated to individually detected sources, 
also contains information about the luminosity function. 
If we assume a luminosity function dN oc Ly~^d\ogL^ 
and a reference luminosity Lr, then the ratio of the con- 
tributions to the total luminosity by sources in the ranges 
(O.l-l)Lr and (l-lO)Lr is given by 
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(1) 



and the ratio of the number of sources in the same ranges 
is 

i?jv(ir) = 10^ = 10i?L(i,.). (2) 

In particular, wc note that for 7 < 0.7 and a luminosity 
function continuing to arbitrarily low luminosities the to- 
tal luminosity of all sources with < is less than the 



luminosity for the sources with (l-lO)Li- In the follow- 
ing we assume that the ratio of counts is proportional 
to the ratio of luminosities, and we use as a reference 
luminosity Lf = 6 x 10^^ erg s~^. 

Following the method described in Johnston & Ver- 
bunt (1996) and Pooley et al. (2002b) we derive from the 
list of detected X-ray sources (Table that 7 = 0.71 
for X-ray luminosities above L^{Q.5 — 2.5) — 4.4 x 
10^9 erg s-^ The K-S probabihty of this value is 93%. 
For K-S probabilities above 10% the slope of the lumi- 
nosity function has 0.47 < 7 < 1.07. 

In addition to the detected sources, we have an ex- 
cess in the core of M4 of about 150 counts, uncorrected 
for absorption, which corresponds to roughly 300 counts 
corrected for absorption in the 0.5-6.0 keV band (which 
is the addition of counts in Xgoft and in Xhard)- In this 
band, the five faintest sources in the core, CX20, CX21, 
CX22, CX26 and CX27 have between 14 and 10 counts. 
A minimum number of sources required to explain the 
excess counts is found by assuming that each source has 
10-14 counts, which gives 22-30 sources. 

The total number of ~ 350 counts from sources in the 
core of M4 with < (five detected sources and the 
excess) is similar to that of the 6 sources in the range (1- 
10) Lr- From our remark following Eqs.^andl^l we see 
that this implies that 7 > 0.7; for lower values of 7 the 
sources with < do not contribute enough counts 
to explain the observed number. On the other hand, for 
7=1 the luminosity function must have a cutoff near 
O.lir, because otherwise the sources with < L^- would 
produce more counts than is observed. With 6 sources in 
the range (1-10) Lr, and for 7 = 0.7 — 1.0 we have 30-60 
core sources in the range (O.l-l)Lr- 

In the core of 47 Tuc there are 14 sources in the range 
(10-100) Lr, which for 7 = 0.7 - 0.8 leads to predicted 
numbers of 70-90 (of which 27 are already detected in- 
dividually) in the range (1-10) Lr, and 350-560 in the 
range (0.1-1) Lr. These numbers are compatible with 
the total number of counts in the core of 47 Tuc from the 
observed sources plus the excess of about 500 counts, 
estimated by Grindlay et al. (2001a). 

In Table |3| we compare the collision numbers and the 
masses of the cores of M4 and 47 Tuc. We expect that 
the numbers of cataclysmic variables scales with the col- 
lision number, and the number of magnetically active 
binaries with the mass. In the range (1-10) Lr the ra- 
tio of numbers of sources in the cores of M4 and 47 Tuc 
suggest that the ratio is set by mass rather than collision 
number. This in turn suggests that magnetically active 
binaries dominate the numbers at these luminosities, in 
agreement with our suggested identifications in M4. 

Comparison with NGC 6397 is also interesting, as its 
collision number is somewhat higher, F = 2.1 (normal- 
ized to the value of M4) while its total mass (as mea- 
sured from its absolute magnitude) is somewhat lower 
than that of M4. Within its half-mass radius, NGC 6397 
contains 9 cataclysmic variables, 4 magnetically active 
binaries and 10 as yet unclassified sources (generally of 
low X-ray luminosity). Within the half mass radius of 
M4 we find 2 or 3 cataclysmic variables, 12 magnetically 
active binaries and 15 unclassified sources, of which an 
estimated 6-12 are background sources. The ratio for the 
numbers in NGC 6397 and M4 of magnetically active bi- 
naries is as expected for a scaling with mass. As already 
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noted by Poolcy et al. (2003), NGC6397 contains rather 
more cataclysmic variables than expected on the basis of 
its collision number. This is reflected in the flat slope 
of its X-ray luminosity function. If the high number of 
bright systems in NGC 6397 is explained by a higher mass 
of this cluster in the past (as suggested by Pooley et al. 
2003), one would expect an accordingly higher number 
of magnetically active binaries, in contrast to what is ob- 
served (assuming that magnetically active binaries have 
a similar evaporation rate as cataclysmic variables, much 
lower than the evaporation rate of single stars). In this 
respect, it is worth noting that the time scale on which 
a binary can be destroyed by a close encounter is short 
in NGC 6397 compared to most other clusters, including 
M4 and 47Tuc (Verbunt 2003, in particular Figure 3b). 
Since magnetically active binaries tend to have longer 
orbital periods then cataclysmic variables, they would 
be easier destroyed by close encounters. This might ex- 



plain the absence of large numbers of such binaries in 
NGC 6397. 
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